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Abstract, de Gennes identified three regimes in the phenomenon of the Langevin 
equation wich includes Coulombic friction. Here we extend and precise this 
phenomenon to a constant external force. 
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1. Introduction 


P-G de Gennes [3] studied the Langevin equation under the influence of a dry friction 
force modelled by the equation 

dv = —-A sgn(v)dt+ VDdB , 

the dry friction force with threshold force A > 0, and D > 0 is the diffusion coefficient. 
Here B is the standard Brownian motion, and sgn{y ) = 1 if v > 0, and sgn{v) = —1 if 
v < 0. Comparing the magnitude of a, A and D de Gennes [3] identified three different 
regimes: viscous, partly stuck and stuck. 

Later Touchette et al. [7] extended de Gennes work by calculating the time- 
dependent propagator of the Langevin equation 


dv 



— a + A sgn{y)\dt + \/~DdB , 


( 1 ) 


which includes a constant external force a G R. 

In this paper, we precise and extend de Gennes’s work to the Langevin equation ([I]) 
and find again the result of Touchette et al. [7] using the trivariate density of Brownian 
motion, its local and occupation times. 
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2. The three different regimes in coulombic friction 


If v(t) is solution of (JTj) , then v(^) satisfies the equation 

dv ————[av — a + Asgn(v)]dt + dB. (2) 

It follows that for large time T the PDF of the velocity v(^) is approximated by the 


stationary PDF 
1 


where 


exp 


Z = — 
2v 


i 

v 2 r 


+ M) 


t- — 2 ?/ 7- — y r + 2y r + y, 

exp(— - )G{ — ) +exp(— -)G(- 


2v yJTV 2v yjrv 

is the partition function i.e. the normalization constant. Here and the sequel 
1 


f + OO 


G(u) = 


V2n . 


. u 2 , D A a 

exp {-—)dv, v=~r, T = — , V = ~- 

2 A a a 


We say that the stochastic process (Vd : D > 0) defined in some probability space 
(0, J P) converges in probability distribution as D —> 0 to the PDF / if for each couple 
l < r of real numbers 

/ r 

f(v)dv, as D —> 0. 

Now we can announce our result. 

1) Stuck regime. If |a| < A, then the velocity v(^) —> 0 as B — )• 0. More precisely 
converges in distribution as D — y 0 to the PDF 
1 — y 2 

—-—exp [-M(l - sgn(v)\y |)]. 

Observe that if the constant force a = 0, then y = 0 and the limit is 
^exp(-|v|). 

2) Partly stuck regime. If |a| = A, then the velocity v(^) —» 0 as D —> 0. More precisely 
we distinguish two cases. 

a) If we consider only the event av () < 0, then 
1 T 

-v (—) 2exp(-2|n|)l [a);< o] as D ->• 0. 


i/ y D' 

b) If we consider only the event av(j^) > 0, then 




2 v 

expf— 


2t' 


[au>0] 


as D —> 0. 


Moreover the probability of the event av(^) > 0 tends to 1 as D —* 0. Hence 
converges to exp(-^)l [a ,, >0 ]. 

3) Viscous regime. If |a| > A then as D —> 0 the velocity v(^) becomes Gaussian 
with the mean (y — sgn(y)r) and the variance vr. More precisely, we have 


*>(£) ~ (y-sgn(y)r) 


v 


1 v 2 
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Observe that the asymptotic mean y — sgn(y)r is the minimizer of the potential 
v + |u| := U(v). 



Figure 1. Three scenarios of the stuck regime with y = 0,0.4,0.9, partly stuck regime 
with y = 1 and r = 1 and viscous regime with y = 3 and r = 1. 

The proof was done in a general case in [3]. For the sake of completeness we recall 
it. It is sufficient to show the case a > 0 i.e. y > 0. 


3. Proof 


3.1. Stuck regime 

2 

We observe that the potential U attains its minimum at v — 0. We have 

' T' 

<r) — 


< 


JTexp(-^)<fo 


V 


/r- P (- 


U(v) 


)dv 


Multiplying the denominator and the nominator by exp(^), and using the change of 
variable - we have 


< 


<|) 

v 


<r) — 


rr 

Ji ex P 

-H(l - sgn{v)l) - 

dv 

r+oo 

J-oo eX P 

[-|u|(i - sgn(v)D - yFg] dv 





















4 


The latter converges to 

fi exp [—|w|(l - sgn(v)^\ dv 


J_T ex P [~\ v \(! - sgn{v )“] dv 
as v —)■ 0, which achieves the proof of the stuck regime. 

3.2. Partly stuck regime 

a) We are going to prove for each l < r < 0 that 

»( 5 ) 


We have 


T r 

< vu ' < r | v (—) < 0) —> / 2exp(2 v)dv. 
v U li 


v(-) T 

P {l < < r \v{—) < 0) = 


v 


D' 


rris 

Jlu eX P 

1 ( 1 (v— t) 2 \ 

A V + It~ )J 

dv 

r° 

J- oc ex P 


I ( V-T ) 2 \ 
” r 2 r > 

dv 


Multiplying the denominator and the nominator by 
variable - we obtain 


ex p(^j) 


, and using the change of 


Ji exp (2v - V"i)) dv 
f- oo ex P ( 2 v -y/”i))dv 


The latter converges to 

jj exp (2 v)dv 


f-oo exp(2v)dv 


as v —y 0, which achieves the proof of the part 1. 
b) We have, for 0 < l < r, 


v (—) T 

P(Z < < r\v(-) > 0) = 

'v D 


rrWu 

Jlyfi ^ 

1 ( . ( v — t ) 2 \ 

l >+ 2r ) 

dv 

f + OO 

Jo ex P 


dv 


ex P( ji;) 


and using the change of 


Multiplying the denominator and the nominator by 
variable ^ we get the proof of the first part of b). 

For the second part we use the same proof and show that P (av(^) > 0) —* 1 as D —* 0. 
3.3. Viscous regime 

The main tool of the proof is the following well known result see e.g. [lj. 

Lemma: Let H be any measurable map such that 

r+oo 


and 


/ exp(— H(v))dv < Too 

-OO 

inf{H(v) : \v — uo| > h} > H(v o) 
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for some Vq and <5 > 0. Then for any 7 > 0, 


/-7 


'[|«-«ol><5] 


—(H(v)-H(v „)) 

V 


dv —> 0 


as ^ ^ 0 . 

Now, let us apply this lemma with H{y) = U(y ) and vq = y — t the minimizer of 
U. We have, for l < r, 

< v{i)-{y-r) < r) = JiS+y~T ex V {~l U ( v )) dv 


XT exp(-it/(n)dr 


We have for v > 0, that 


U{v) — U(y — t) = 


(v-(y- t)Y 


2 r 


If l > — 00 , then for small 17 we have 


rry/u+y—T 

r 1 1 

dv 

rr^/u+y—T 

[ 1 nl 

/ exp 

/ ly/u+y—T 

— (t/(u) - U(y - r)) 

V 

V^ = 

/ exp 

J ls/v+y-T 

2 rzy ( v ( y r )) 


dv 


r v 2 

= Ji exp(——)dv. 


and then 


< < r) = r exp( _^ ) * . 

1 G ~ ’ Ji K 2r y/2nr 


If l = — 00 , then 


PTy/u+y—T 

1 

/ exp 

/ — 00 

-~(U(v) - U(y - r)) 


dv 

-^ = (l) + (2), 


where 


( 1 ) = / exp 
J [w<0] 


v 


dv 


—y 0 , 


( 2 ) = / exp 

J [ 0 <v<r^/L'+y— t] 


~-(U(v) - U(y - t)) 


dv 


From Lemma (13.3ft the term (1) converges to 0. By the change of variable z = u ^ T ' ) . 
the term ( 2 ) 


( 2 ) = 


l[-^<z< r] 


eXp ( _ 2 r) dZ 


converges to exp {—^)dz. By taking r = + 00 , we get 


^+00 


exp 


— (U(v) - U{y - r)) 


dv f + °° z 2 

, ->■ / ex P {~^~)dz, 

'v 2 r 


and then 


< 


«(S)-(y-T) 


r v 2 , dv 

* r W, 


which achieves the proof. 
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4. Time-dependent propagator 

Now we drop the coefficient \ in (JT]) and we discuss the calculation of the time-dependent 
propagator of 

dv = ~[av + a + A sgn{y)\dt + \/~DdB. 

Using the equality of the laws or the probability distributions of (y/DB(B) and ( B{t )), 
we derive that 

Law(v a,a ' A ’ D (t)) = Law(v^'^'Ti )l (Dt)). 

Hence the propagators p a ’ a ’ A ’ D ( v ^ f | Vq ^ q), 1 1 VQi o) respectively of v a,a,A ’ D (t) 

and satisfy the relation 

p a,a ’ A ’ D ( Vj t | v 0 , 0) = Dt | v 0 ,0 ). 

Hence, it sufficient to study the case D = 1. 


5. Time-dependent propagator for a = a = 0 using local occupation time 


We denote by P and P „ 0 the probability distribution respectively of the trajectories 
s G [0, t] —> v(s) of the solution of (JT]) and the Brownian motion starting from Vq. 
Under the probability distribution 

tA 2 ' 


exp 


sgn(B s )dB s 


d^vo ■= fsgn(B)dF. 


VQ 


the process ( B(s ) : s G [ 0 ,t]) is solution of the equation 
dv = —A sgn(y)dt + dB , u(0) = v 0 . 

We simplify the stochastic integral f* sgn(B s )dB s using Tanaka formula [ 6 ] 
\B t \ = |n 0 | + [ sgn(B s )dB s + 2 L t . 


Here the local time 


L t = lim — 

£->o 4e 


L [ffia|<£] 


ds 


1 

2 


It follows that 


5(B s )ds. 


sgn(B s )dB s = |n 0 | - \B t \ + 2 L t . 


Now, 


f S gn{B ) = exp ( A(|n 0 | - \B t \ + 2 L t ) - 


tA 2 


The densities of v[t) and the Brownian motion B{t) are related by 


p(v,t | v 0 ) = E. 


'^0 


5(B t - v ) exp(A(|n 0 | - \B t \ + 2 L t ) 


tA 2 , 


(3) 


The latter formula is also known as path integral representation [2]. Hence the law of 
the solution v(t) is given by the law of ( B tl L t ). 
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5.1. Density of Brownian motion and its local time 


Set r t = Jq 1 [B s > 0]ds, and 

I I 2 

\v\ V 

h(s,v) = -/==^= ex P(-7r)’ s > 0, n G M. 

V ZS^Tl 

Karatzas and Shreve |5j have calculated the probability density P V0 (B t G db, L t G dl,T t G 
dr) := pt(dv, dl, dr | no) of (B t , L t ,Y t ) as follows. For no > 0 we have 

P Vo (B t G db, L t G dl, T f G dr) = 2 h(r, l + vo)h(t — r,l — b)dbdl, b < 0, 
P„ 0 (£?t G db, L t G dl, r 4 G dr) = 2h(t — r, l)h(r, l + b + no )dbdl, b > 0, 
P„ 0 (Sj G db,L t = 0 ,r t = t) = cu(v 0 ,b,t), b> 0,n o >0, 


where 


u}(v 0 ,b,t) = 7 t (b - vq) - 7 t(& + n 0 ), 


7t(w) 


1 r ^ 

. expl- 

1 2 t ; 


We derive the joint distribution of ( B t , L t ) under P „ 0 with no > 0: 

(21 + n 0 - b) 2 


P V0 (B t G db, L t G dl) = 2 ( i—^==^~ exp 

P V0 (B t G db, L t G dl) = exp 

\Z2nt 3 


2 1 

(21 + n 0 + b) 2 
2t 


dbdl, b < 0 , l > 0 , 


dbdl + u(vo,b,t)5(l), b>0,l>0. 


p(v,t\v 0 ) 


Now, we calculate the density of the solution (CTJ) as follows. If n < 0, then 

( tA 2 ' 

5(B(t) - v) exp ( A(n 0 - \B t \ + 2 L t ) -- 


= E. 


'vo 


:= exp |^A(n 0 + n) 

= 2 exp I A(n 0 + n) — 


tA 2 

2 

tA 2 


= exp(A(n 0 + n) 


tA 2 


E„ 0 [S(B(t) - n)exp(2A L t )} 

" + °° ( 2 / -v+_ ,„) _ (2l-v + v 

, V2f% 2 1 

+ ” (i - t, + Vexp(Ai- ( ^A 1, ° )2 M 


V2D 


7T 


21 


tA 


2 r 


= exp(A(n 0 + n)-—) 


-y/2tTT 


exp(- 


(n 0 - vf 


"+00 


2 1 


+ A exp (A/ - 


(l - n + n 0 ) ; 
2 1 


dl 


■\Z2tn 


1 / tA 2 . . . .. . (n 0 — n ) 2 

exp(-—) exp(A(n 0 + n)) exp(- v u y ) + 


\J 2tn 

tA 2 

A exp(A(n 0 + n)-—) 


r*+oo 


exp(AZ — 


2t 

(1 + vq-v) 2 . 


dl 


2t y/2tn 

After some calculation we obtain 

f + °° /a, (/+ no — n) 2 dl A 2 t / A/ .. ,v — no + At, 

/ exp(AZ---) = exp( —-) exp(A(n — vq))F(- 


'0 


2 1 


y/2tir 


Vi 


where F(v) — f_ 


v exp{-\-) 


—du. Finally for no > 0, n < 0, we have 
tA 2 , , , , . , ,, „,n — no + At, 


V2it 


p(v, 1 1 n 0 ) = exp(——) 7 t(n 0 - v) exp(A(n 0 - v)) + F(- 


Vi 


Aexp( 2 An). 















































If v > 0, then 
p(v,t\v 0 ) = E. 


'’0 


S(B(t) - v ) exp (A(n 0 -\B t \ + 2L t ) - 


tA 2 . 


:= exp (A(v 0 - v ) 


tA 2 


2 
t 

= exp (A(n 0 - n)-—) 


)E„ 0 5(B(t) - v) exp(2AL t ) 


f+OO 


(21 + n + no) 


V2h 


7T 


exp(2A/ — 


(2/ + v + no) 2 
2t 


)dl + u(v 0 ,v,t) 


= exp(A(n 0 - n) 
+ exp(A(n 0 - n) • 


tA 2 . 

9 ' 

tA 2 , 


u(y 0 ,v,t) 
r°° (i+v+v o) 


^2f3 


7T 


( a ; + V + n 0 ) 2 

exp (A/-—- )dl 


... . tA 2 tA\r 1 , (u 0 + u) 2 , 

= exp(A(n 0 - n)-—)u;(u 0 , n, t) + exp(A(n 0 - n)-—) - exp(- 


2 Ay/2t^ 


2t 


^+00 


+ A 


exp (A l 


(l + v + n 0 ) 2 dl 


2 1 


\j2tix- 


... . tA 2 . ,1 tA 2 , . (u 0 + n) 2 , 

= exp(A(n 0 - n) - — )cn(n 0 , n, t) + -^= exp(-—) exp(A(n 0 - n)) exp(-—-, 


+ Aexp(A(n 0 - n) - -^-) 

From some calculation we obtain 
, ‘ +00 (/ + n 0 + n) 2 . dl 


n+0 ° (/ + n 0 + n) 2 N dl 

exp(A l - ' 


/“too 

/ exp(A/ - 

Jo 


A 2 t _ 

2t tfe? = exp( —> eXp( “ A( '' + B ° ))F( ~ 


2t y/2t~K 

At - (n + n 0 ), 


Finally if n > 0, n 0 > 0, then 

, | . , , . tA 2 . , 1 . tA 2 (no + n) 2 . 

P(v, 1 1 n 0 ) = exp(A(n 0 - n)-—Mn 0 , n, t) + -^== exp(A(n 0 - n)-—) exp(-—-) 


+ Aexp(-2An)F( 
1 


At - (n + n 0 ). 


Vt ' 

... x tA 2 . (vo — n) 2 . , „ . . At — (n + n 0 ). 

exp(A(n 0 - n) - —) exp(- v 0 ; ) + A exp(-2An)F(--A) 


■\Z2tn 2 21 y/t 

= ^exp(A(n 0 + n) - ^-)y t (v - n 0 ) + F( — — ^_ + ^ )^ Aexp(-2An). 

Finally we have for u,no G M, that 

p(v, 11 n 0 ) = g(n, 11 n 0 ) exp(-2A|n|) 

where 

g(n, 11 n 0 ) = A ^exp(A(|n 0 | + |n|) - ~ v 0 ) + F(— — 

Observe that q(v,t | no) is symmetric, i.e. q(v,t | no) = q A (vo,t\v). In the language of 
linear diffusion m(v) = exp(—2A|n|) is the speed measure of the linear diffusion ([3]). 
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6. The case a ^ 0 


In this case the probability distribution IP of the solution 
dv = —[A sgn(y) + a]dt + dB, u(0) = vq, 

is also absolutely continuous with respect to P„ 0 (the probability distribution of the 
Brownian motion starting from uq). We have 


dP 

dP v 


( B ) = exp(- / ( Asgn(B s ) + a)dB s - - / (Asgn(B s ) + a) 2 ds). 


' o 


v 0 J 0 

After some calculation we have 

- / (A sgn(B s ) + a)dB s = A(|u 0 | - \B t \ + a(u 0 - B t )) + 2A L t , 
Jo 

j (A sgn(B s ) + a) 2 ds = (A 2 + a 2 )t + 2aA(2T t — t). 


It follows that 


p(v,t | vq) = exp 


A(|u 0 | - |u| + a(v 0 - v)) - 


(A — a) 2 t 


Ed 0 [5(B t - v ) exp(2ALj - 2aAT t )] . 


Then p(v,t | u 0 ) is calculated using the trivariate probability distribution p t (db,dl,dr) 
of (B t ,L t ,T t ) as follows: 


p(v,t\v 0 ) = exp 


A(|u 0 | - |v| + a(v 0 - v)) - 


(A — a) 2 t 


r»+CXD ft 


exp(2A(/ — ar))p t (v, dl , dr). 


>o J o 


7. The general case 

Similarly as above the density of the solution of 

dv = —[av + A sgn(v) + a]dt + dB, u(0) = v 0 , 

is 

p(v,t | vq) = exp 


a /i i ii / \s (A — a) 2 t at 

A (ko| - M + a(vo ~ v)) ---+ — 


+00 ft f~\~ oo /»+OO f~\~ oo 



a 


exp (2A l — 2aAr — —b 2 — qA|6i| — aab\)p t (v, dl, dr, db\, d\b\\, db 2 ), 


'0 J 0 JO JO J—oo 

where pt(db, dl, dr, db\, d\b\\, db 2 ) is the probability density of 

(B t ,L t ,T t , j B s ds, I \B s \ds, I B 2 ds). 

Jo Jo Jo 


8. Conclusion 


We have precised and extended the three different regimes of the Langevin equation 
which includes a viscous friction force, a Coulombic friction and a constant external 
force. Moreover we find again its time-dependent propagator using the density of 
Brownian motion, its local and occupation times. 
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